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Abstract 

We shall present a new strategy for handling mean field limits of quan- 
tum mechanical systems. The new method is simple and effective, ft is 
simple, because it translates the idea behind the mean field description of 
a many particle quantum system directly into a mathematical algorithm. 
It is effective because the strategy yields with lesser effort better results 
than previously achieved. As an instructional example we treat a simple 
model for the time dependent Hartree equation which we derive under 
more general conditions than what has been considered so far. Other 
mean field scalings leading e.g. to the Gross-Pitaevskii equation can also 
be treated [151 PH] . 

1 Introduction 

The dynamics of a quantum mechanical many body systems with interaction 
can sometimes be well approximated by an effective description in which each 
particle moves in the mean field generated by all other particles. Derivations of 
such mean field equations from the microscopic N body Schrodinger evolution 
are usually done for the reduced one particle density and are naturally based 
on hierarchies [M 13 1 SI UJ • 

In the recent years alternative methods have been succesfully used to derive 
the Hartree equation from microscopic dynamics. One approach was developed 
by Frohlich et al. using dispersive estimates and counting of Feynman graphs 
[7] . Another one was introduced by Rodnianski and Schlein jTS] ■ They focus on 
the dynamics of coherent states, inspired by a semiclassical argument of Hepp 

We present here a new method for deriving mean field descriptions which 
is in particular when deriving the Gross-Pitaevskii equation simpler and more 
effective as it yields more general results with greater ease. For concreteness of 
the presentation we consider a simple model leading to the Hartree equation. 
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We consider a Bose condensate of N interacting particles when the external trap 
described by an external potential A* — is time varying, it can for example be 
removed. We are interested in solutions of the iV-particle Schrodinger equation 
in units H = m = 1 

i¥ N = H^ N (1) 

with symmetric initial wave function "I^ when N gets large and the interaction 
gets weak with N: The Hamiltonian 

TV N 
j=l l<j<k<N j=l 

acts on the Hilbert space L 2 (M. 3N ), and (3 £ R determines the scaling behavior 
of the interaction. Usually v^f scales with the particle number such that the 
total interaction energy scales in the same way as the total kinetic energy of 
the N particles. This means that the L 1 -norm of is proportional to iW 1 , 
for example v^(x) = N^ 1+3 ^v(N^x) for a compactly supported, spherically 
symmetric, positive potential v € L°° . Thus the total interaction energy is for 
sufficiently smooth wave functions ^ of order N. For positive (3 the support 
of the potential shrinks with N. As long as (3 < 1/3 the interaction potentials 
overlap and the mean field approximation is heuristically clear. 

For 1/3 < (3 < 1 the interactions get more 5-like and do not overlap. But in 
this case the wave function \I> develops on the scale of the potential a structure 
around the centers of the interactions to keep the energy low. If the energy 
of \E' is controllable the interaction effectively still behaves like a smeared out 
interaction with moderate scaling behavior and the mean field argument still 
holds [14]. 

The trap potential A 1 does not depend on N. Hn conserves symmetry, i.e. 
any symmetric function evolves into a symmetric function ^> l N . 

Assume for the sake of simplicity for now that the initial wave functions 
is a product state = JljLi <P°( x j) where tp° e L? . 

In the mean field limit the product structure survives during the time evo- 
lution, i.e. for N large 11,3=1 l P t { x i)i but every particle moves in the 
average field of all other particles (mean field) so that ip* solves the Schrodinger 
equation 

=(-A + A* + (3) 

( with ip° as above). The "mean field" V t depends on ip l itself, so ([3]) is a 
non-linear equation. 

Our new strategy revolves around the meaning of ~ IljLi ^'(^j) f° r N 
large: "Most" particles behave in a good way and the term on the right has 
"mostly" product structure while only "few" particles will behave badly and 
will become entangled. We shall therefore introduce a biased counting of good 
and bad particles yielding a counting measure a(t) such that for a ~ most 
particles are good. The algorithm will then produce an equation for a(t) which 
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shows that if a(0) s» then a(t) 0. This result is then easily generalized to 
non product initial states and also to an assertion about the reduced one particle 
density matrix, which is usually the way limits are phrased in [5J [5J 1151 116] : 

fh N ( x , V) '■= J *at( 2; > • • • ) x N )^ N (y, x 2 ,..., x N )d 3 x 2 ■ . ■ d 3 x N . 

[i 1 N converges to \tp t )(tp t \ in trace norm. Such results are usually based on a 
hierarchical method analogous to BBGKY hierarchies. 

A warning: One may be inclined to think of s» Y[j=i V^&j) m L 2 sense. 
That is however false: Assume that 

N-l 

= \\ (f(xj)(f L (xjf) f or some (p^-Lip ■ (4) 
i=i 

Then of course ^at-L YljLi <f( x j)> so m £ 2 -sense ^jy is far away from J\ip. 
Nevertheless clearly "most" particles of ^jy are in the state tp. That is what we 
must focus on. 



2 The strategy 

We wish to control the number of bad particles in the condensate. The idea is 
simply to count (but in a biased way) the relative number of particles not in 
the state tp*. This then leads to the following definition: 

Definition 2.1. For any tp 6 L 2 (R 3 ) and any 1 < j < N the projectors Pj ° : 
L 2 (R 3N ) -> L 2 (R 3N ) and qf : L 2 {R 3N ) -> L 2 {R 3N ) are given by 

P J* N = tp( Xj ) J p*(x J )^ N (x 1 ,...,x N )d 3 x ] V* N tL 2 (R 3N ) 

and qj = 1 — pj. We shall also use the bra-ket notation pj — \tp(xj))(tp(xj)\. 
Furthermore we define on L 2 (R 3N ) for any < k < N the projector 

/ k N 

P V (ll-/: n *? # 

sym 

The index sym means the following: For any < k < j < N consider the set 

N 

A k ■■= {(ai,a 2 ,...,a N ) : aj £ {0,1} ; J~] a 3 - = k} . 

i=i 

Then 

N 

, 1 — a-i 
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With k at hand we can define an object which "counts the number of 
particles which are not in the state tp in . 

Definition 2.2. Let ((•,•)) be the scalar product on L 2 (M. 3N ). We define for 
any N £ N and any function n(k) : {0, 1, . . . , N} — > the functional '■ 
L 2 {R 3N ) x L 2 (R 3 ) as 

N 
k=0 

Note that uf =0 A = {0, 1} N , hence Y,k=o P N,k = 1 - 

Let us explain a bit more what c<n does. Choose n(k) — k/N. Then 
that part of where k of the N particles are not in the state (i.e. 
{^n, fc^w))) is given the weight k/N. Hence ajv "counts the relative number 
of particles which are not in the state ip*" . 

The important role ajv plays lies however in the fact that the function 
n{k) can be chosen more appropriately depending on the particular problem 
(see section 12.41 below) . The choice will be dictated by the requirement that 
ctN^N, implies convergence of /^*™ to in trace norm. 

2.1 Control of a 

Let us use the shorthand notation a(t) := aN(^%, Our goal is to prove 
that limjv->oo ct(t) = if limj\r->oo ct(0) = 0. Using Gr0nwall it is sufficient to 
show that 

\a(t)\ <Ca{t)+o(l) . 

Now 

N 

a N (K, f") ■= [Hn K f , n(k)Pl k ]^ N )) (5) 

k=0 

where in view of © H™ f is the iV-body mean field Hamiltonian H N = 

Note that many of the terms appearing in (c.f. ©) and cancel in 
the difference Hn — H 1 ^ . All that remains is the interaction potential minus 
the mean field potential. 

It is important to note that in this algorithm no propagation estimates on 
^f f N are needed. It is possible to estimate the right hand side of ([5]) 

N 
k=0 

in terms of a(^N, y), o(l) and the energy of uniformly in and if. 
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2.2 Advantages of the method 

• No hierarchies appear 

• No propagation estimates are needed 

• The freedom of choice of the weight n(k) gives flexibility in the proof (see 
section 12.41 below) . 

2.3 Convergence of the Reduced Density Matrix 

An important feature of the functionals ajy defined above is that 
limAT^oo ajv(^iV) <p) =0 implies convergence of the reduced one particle density 
matrix for many different weights n(k) and vice versus. In particular we can 
state both our condition (a(0) — > 0) and our result {a(t) — > 0) on the respective 
convergence of the reduced one particle density matrix instead. 

Before we prove this equivalence note first that the weights n(k) = k/N has 
the special property that 

N b N 

^-E^^E^r w 

k=0 j = l 

To see this we use that ^Zk=o-^ > Nk ~ ■"•> hence it suffices to show that for 
any < k < N 

N . N 
3=0 3=1 

The left side equals £p* fc . Recall that P* k = {]\ k J=l qJUjLk+iPkUm- Mul- 
tiplying this with Xw=i if yields a factor k. Thus also the right hand side 
equals ^P% >k . 

Because of ^ among the different weights n(k) — k/N is easiest to han- 
dle. Therefore we shall show first that \imN^ o(('^N,n lp ' ( &N} = is equivalent 
to convergence of the reduced density matrix and generalize to other weights 
thereafter. 

Note that convergence of fif N to \<p){<p\ in tr ace norm is equivalent to con- 
vergence in operator norm, since |<p)(<p| is a rank one projection |15j . Other 
equivalent definitions of asymptotic 100% condensation can be found in |12j . 

Lemma 2.3. Let j > and n(k) = k/N . Let <p G L 2 and let ty N € L 2 

e sy 

(a) 



be symmetric. Let /if N be the reduced one particle density matrix o/^jv- Then 



(b) 



lim ((^jv, ^^tv)) = ^ l™ Mi N = I^K^I ^ n operator norm. 

N— >oc N— >oo 



lim {9 N ,n v ^ N }) =0 <S=> lim {V N , (n^f V N )) = . 

N—tOO N— fOO 
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Proof. Note first that using © and symmetry of 

{9 = \\q(y N \\ 2 ■ 

(a) "^>" Let 



lim ({^ N ,n^ N } = 0, 

jV— s-oo 

i.e. limAr^oo ||gf $jv|| = and lim7v_ J . 00 ||p^jv|| = 1. Note that 
/if N (x, y) = 1 $ N (x,x 2 , ■ ■ ■ , x N )^(y, x 2 , ■ ■ ■ , x N )d 3N ~ 3 x 



l pf , t> N (x i X2, ■ . . 


, x N)Pi^* N (y,x 2 , . . . 


,x N )d 3N ~ 3 


f J qf^N(x,x 2 , ■ 




..,x N )d 3N 


+■ J Pi^n{x,x 2i ■ 


■■, x N)qf^* N (y, x 2,- 


..,x N )d 3N ' 


■f J qf 1 $'N(x,X2, ■ 


■ ■>x N )qf^* N (y,X 2) ■ 


..,x N )d 3N ' 



The hrst summand equals HpfWjvll 2 \(f)(<p\- The second and third have op- 
erator norm \\qf ^>n\\ Wpf^wW an d the fourth has operator norm H^f'I'ivll 2 
and hence go to zero. 

(a) "-4=" Assume that /i*" — > in operator norm as N — > oo. It follows 

that lim A r^ 00 ((/3, /if N tp) — 1. Writing /if N like in J7]) and using that 
qftp(xi) = the second, third and fourth summand are zero. Therefore 

lim \\pfV N \\ 2 = l. 

N— foo 

Using now pf + qf = 1 it follows that limAr^oo Hgf^jvll = 0. 
(b) For (b) we show that 

lim (* N , (ny * N } = lim {9 N , (r^f * N } = (8) 

N— >oo N— >oc 

for any j, I > 0, which is equivalent to (b). 

Let limjv-foo^Ar, (n v y^N} = for some j > 0. We shall use the abbre- 



viation 



N / j \ j 



fe=0 



Let fcjv be the smallest integer such that (%■) < v^iv. It follows that 
(fc*±l) j > ^ an d thus ll^fc^Jvll 2 < V^- Hence 

N / h\ l kN / h\ l N 

El -I) n^-ii 2 * E l iiv^+EW 



N i m jv,* - ^ y iv . 

fe=0 v 7 fc=0 v ' k N + l 



Nil" 



< 



G 



Thus limjv-voo^JV, {n^f *jv}> = and © follows implying (b). 

□ 

Remark. Similarly one can proof that limAT-j-oo^iV) (n^Y'&N} = for j G M + 
implies convergence of the reduced k-particle density matrix for any fixed k < oo. 
Also note that for any m{k) < n(k) 

N N 

lim «^,£m(fc)P£ fc *jv)) = lim «*jv, £ = . 

iv— >-oo * — ' iv— ^oo * — ' 

fe=0 fc=0 

From this follows that limjv-^oo ctjv( v l'jv, <p) = implies convergence of the re- 
duced one particle density matrix for any weight dominated by {k/N) 3 for some 
positive j . 



2.4 The role of the weight 

An important feature of this new method is the freedom of choice for the weight 
n(k). In the instructional example below we will explicitely deal with the sim- 
plest scaling behaviour which is f3 = 0. In this case we get a Gronwall-type 
estimate for a(^>t, 'fit) for many different weights in particular for any 

weight that can be written as n(k) = (k/Ny for some positive j. We shall 
choose n(k) — k/N below, which is due to ([6]) the most convenient choice. 

But for other situations other choices for the weight are more appropriate: 
When dealing with scalings (3 > 1 one can either 

(a) use high purity of the condensate to control -j^a^t, 'fit) although the 
interaction gets 5-like. 

(b) control some of the kinetic energy and use smoothness of ^. 

Both ideas can be worked out, and fundamentally depend on the choice for the 
weight. Let us explain 

(a) is worked out in |14j where we have to restrict ourselves to scalings /3 < 
1/6. Still the result is interesting since it is so far the only derivation of the 
Gross Pitaevskii equation without positivity condition on the interaction. 
Here we use scalings of the form 



n(k) 



k/Ni iik<N"i 
else 



with < 7 < 1. Depending on 7, a — ¥ stands for a different purity 
of the condensate: For 7 = a — > is equivalent to L 2 -convergence of 



^{xi, . . . , Xn) against the full product Yij=i f( x j)> f° r 7 = 1 we have the 
"old" weight n(k) = k/N. For "large" /3 one needs high purity of the 
condensate to derive the Gross-Pitaevskii equation: For f3 close to 1/6 7 
has to be chosen close to 0. 
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(b) is worked out in [T3] and gives — assuming positivity of the interaction — 
good results for all scalings < /3 < 1: Doing the estimates one arrives for 
any weight n(k) — {k/N)i with positive j roughly at the following formula 

j t a(^ t ,ip t ) < C.ia^t^t) + o(l) + \\V iq ^ t \\ 2 ) . 

Now ||Vi<7jf 'Pill can be controlled using conservation of energy and split- 
ting up (y>t,H mf ^ t )- It turns out that for weights n(k) = (k/N) 3 with 
j < 1/2 one can show that || Vigf ^t\\ 2 < Ca(^>t, ft)- Choosing for exam- 
ple the weight n(k) — y/k/N one arrives at a Gronwall type estimate for 
a(^ t ,ip t ). 

A similar idea can be used in the Hartree-case (i.e. [3 = 0) when dealing 
with interactions with strong singularities (see [TU]). 

3 A simple example: The Hartree equation 

For the rest of this paper we will apply the algorithm to the case [3 = for ease 
of presentation. Nevertheless our result here is already better than previous 
results [2] |51 [IS] in the sense that we can handle stronger singularities than 
Coulomb in the interaction. 

With more technical effort it is possible to treat, even stronger singularities 
in the interaction [10]. Other scalings ((3 > 0) are treated in [13] and [14] . 

For the scaling [3 = the mean field potential is simply the convolution 
v*]^] 2 , i.e. ([3]) becomes the Hartree equation 

idY =h H ip t =: (- A + A'+t^l^py . (9) 

We choose for the weight n(k) = k/N. Assuming that aN(^%, <f°) — > as 
N — > oo we show that a/v^jv; (p*) — > 0. 

Theorem 3.1. Let v G L 2r for some r > 1. Let A 1 be a time dependent poten- 
tial. Assume that for any JV S N there exists a solution 'J'^y of the Schrddunger 
equation (QP and a solution tp* of the Hartree equation with \\(p \\2 S — o° for 
s = -t—r . Then 

r—l 

MK, <P*) < ^ ^(n, <P°) + (e f o CTdT 1)N^ , 
where C* := 10| | w|| 2r- 1 1 II 2s • 

Remark. (a) //Halloo < oo for all t < oo we can handle interactions v € L 2 . 
The Theorem generalizes the results in the literature to time dependent 
external fields. 

(b) Assuming limTv-^oo ajv^jv, l fi°) ~ an d Jo°° II 2s<^t < oo the Theorem 
implies that liniAr^oo aiN^^, ip l ) = uniform in t < oo. 
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(c) There is a lot of literature on solutions of nonlinear Schrddinger equations 
(see for example ]3j) showing that our assumptions on the solutions of the 
Hartree equation can be satisfied for many different setups. 

As mentioned above we do not need any propagation estimates on ^ l N . To 
emphasize this we prove a stronger statement than the one in the Theorem. We 
define the functional j N : L 2 (R 3 ) <g> L 2 {R 3N ) -> K by 

7w (*at, ip) := -i({$ N , [H N - H§, n^ N ) (10) 

where 

ff^:=^-A i +A t (^) + («*|^| 2 )fe) 

3=1 

is the sum of Hartree Hamiltonians for each particle. Recall the ip 1 solves ©, 
so aw and are such that 

a N (* t N ,<p t )= lN (y t N ,ip t ) . (11) 

77v(vt , Ar, f) will now be estimated in terms of ct]\r(^N, f) uniformly in \1/jv 
and ip. 

Lemma 3.2. Let v G L 2r for some r > 1. Then 

\7n{Vn, <p)\ < 10C V (a N ($ N , if) + N- 1 ) 
for all ^jv € L 2 and all if € L 2s with s — -^—^and C v := || f|| 2t- || vll 2s - 

From this Lemma the Theorem follows in view of (| 1 1 [) using Gr0nwall: Writ- 
ing C" := 10|M| 2 ,||^|| 2s and / := a N (^ N , f*) + N' 1 we have |/*| < C*/*. It 
follows that /* < e^o CTdT /°, i.e. 

ajv(^, y*) + N- 1 < eti ° TdT (a N (^ N , <p*) + A^ 1 ) 
which proves the Theorem. 
Proof of Lemma \3.2[ Recall that 

N 

H N - H% = ^2 v n{xh ~ xi) - ^2(v* \f\ 2 ){xi) 

l<j<k<N 1=1 
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and n^j = N 1 Ylf=i if - Using symmetry of and 1 — pf + qf it follows 

JV N 

Jn(Vn,<p) = -iN- 1 ^{^ N ,^v N {x k -x l )-^£ j v^\if,\ 2 {x l ),qf\^ N ) l 

j=l Kk 1=1 

N 

= -i((y N ,[^2v N (x k - an) - (v * \ip\ 2 ) (an), gjf ]*at» 

fe=2 

= -i|* JV ,((JV-l)^(a;2-a;i)-(t;*M 2 )(a;i))gJ'*^) 
+i«*jv,gf ((JV - - an) - («* M 2 ) (an))*iv» 

= -*|* JV ,gf((JV-l)i;^(s2-si)-(i;*M 2 )(a;i))g]['*^) 
-i«*iv,Pr((JV - l)v N (x 2 - xx) - («* M 2 ) (aJi))gST*jv) 
+i|* JV ,gJ' ((JV - l)ujv(ar 2 - an) - (v* \tp\ 2 ) ( Xl ))q^ N } 
+z((*at, qf({N - 1)ujv(^ ~ an) - («* M 2 ) (xO^jv)) 

= 23(((* JV , P n(^- 1 W(z2-2:i)- («*M 2 )(si)K*JV» . 

With 1 = p*2 + q 2 we have that 

= 2Z((* N ,p^((N-l)v N (x 2 ~x 1 )-(v*\ l p\ 2 ))( Xl )q?p^ N ))) (12) 

+23 - i) VJV (x 2 - Xl ) - (v* \<p\ 2 ) ( Xl )) q r q ^ N ))) 

+23 (((H> N ,pUi((N - l)v N (x 2 - Xl ) -{v* M 2 ) (xi))qfp^ N ))) 
+23 (i*N,pU$((N - 1)M*2 - *i) " («* M') M«*iV>>) ■ 

Since for any selfadjoint v4 on L 2 (M. 3N ) the operator Pi<7 2 -^QiP^ ^ s invariant 
under adjunction with simultaneous exchange of the variables x\ and x% , we see 
that the third summand is zero. 

Let under slight abuse of notation for any j > 

N ( k \ - j 



k=l 



It follows that (n^) J (n v )~ J +P* Q = 1, thus for any 1 < j < N 

(nvy(n*)-jqf = qj . (13) 

Defining 

V(xi,x 2 ) ■= (N - l)v N (x 2 - Xl ) - (v* \(p\ 2 ) (on) 
and using (| 13[) on the second summand of (fT2|) we get 



|7jv(*jv,v)I < 2\(* N ,p*p%V(x 1 ,x 2 )q?p 2 °* N }\ (14) 
+2|((vI/ Ar! ^y( a;i!a;2 )(^) 1 / 2 (^)- 1 /2^^v|/, v ))| 

+2|((*^,^g| I y(a ;i ,a ;2 )grg 2 ^ J v»| . 
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The first summand is the most important. It becomes small because the 
interaction is well approximated by the mean field potential. Recalling the 
notation p^ = \(p(x2)){<p(x2)\ and the scaling behavior of vn — N~ x v it follows 
that 

p$V{x 1 ,x 2 )p2 =(N - l)p%v N (x 2 -x x )p% -p% (vk\(p\ 2 ) (xi)p% 

= (1 - N-^ix^iipixzMxz - Xl )\<p(x2))(<p(x2)\ -pi (v * M 2 ) { Xl ) 

= (i - n- 1 )]^)) (« * m 2 ) {x x ){<p{x2)\ -fi{v* M 2 ) On) 

= -N~ l pt (v*\<p\ 2 ) (xx). 
Hence the first summand in (|14p equals 

< 27V- 1 !! («*M 2 ) ( Xl ) P fp^ N \\ \\ q f P ^ N \\ 

< 2N- 1 \\(v*\i P \ 2 )(x 1 )p^ N \\ 

= 27V- 1 (((* J v,|^ 1 ))(^(x 1 )|(^|^| 2 ) 2 ( a ; 1 )|^(x 1 ))Mx 1 )|vI/ Ar ))) 1/2 (15) 
With Young's inequality 

II ( V *| V | 2 ) 2 H,. = Hv^M 2 !!^ < (||v 2 lli]]«||*.) 2 = Hvlllr - 

Using Holder inequality recalling that 7 + 7 = 1 (fT5|) is bounded by 

2^- 1 (||(^|^| 2 ) 2 || r ||^ 2 || s ) 1 / 2 <27V- 1 (|| U || 2 r ||^|| 2 s ) 1/2 =27V- 1 ^. (16) 

Next we estimate the second summand of (fT4|) . Using ||<^|| = 1 and the 
scaling vn = iV _1 u we get 

Sup ||V(.,a!2)||r <(N~ l)\\v N \\ r + ||«*M 2 ||r < ||»||r + \Hr\\<pf = 2||«|| r . 

Thus we have in operator norm 

\\pfV 2 (x 1 ,X2)pf\\o P < WfX sup ||U 2 (-,x 2 )|| r < 4||^||I s || U ||l r < 4(C") 2 • 

(17) 

Going back to (| X4[) and using Schwarz inequality the second summand there is 
bounded by 

2\\(n'0) 1/2 V(x 1 ,X2)pM*N\\ \\(nn- 1/2 qfq^N\\ ■ (18) 



11 



Using symmetry 

\\(nn 1/2 V(x u x 2 )pfp^ N f 

N 

= N-iY,i*N,pf P2 °V(x u x 2 )qfV(x 1 ,x 2 )p* P $*N} 

3=1 

= ^j^(* N ,p'?p%V(xi,x 2 )q$V(x u x 2 )pfp 2 >* N )) 

+^(* N) p$p 2 °V(x u x 2 )qfV(x 1 ,x 2 )pfp 2 °* N )) 
< H9 N f \\piV 2 {x u x 2 )pf\\ op +^\\pfV\x u x 2 )p^\\ op 

<4(Cn 2 a N (* N ^) + ^p^ . (19) 
Using symmetry and © 

N(N-l)\\(nn- 1/2 qfq^N\\ 2 = N(N - 1)(* N , (n*)" 

AT 

< ^«* JV! (^)- 1 g J^ JV » 

where we used (fTF)) to get the last line. This and (fT9|) yield that (fT8| (i.e. the 
second summand in (fT4|) ) is for AT > 1 bounded by 

3C ,v ^/^(*^)(a J v(*iV, ¥>) + N- 1 ) 1 ' 2 < 6C*{a N (V N , if) + N^ 1 ) . (20) 
Using Schwarz inequality and ()17|) the third summand in (|14[) is bounded by 

2||U(* 1 ,x 2 K^vM II^^H < 2(\\pfV 2 (x 1 ,x 2 )pf\\ op ) 1/2 \\q$* N \\ 2 

< iC^a N (^ N , <p) . (21) 



The bounds ([16]). (|20|) and (|2Tj) of the three summands in ((T4J) imply the 
Lemma. 

□ 
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